function [x_noise,std_dev] = addNoise( x, level )
% function [x_noise, std_dev] = addNoise( x, level )

std_dev = level*( max( x(:) ) - min( x(:) ) );
x_noise = x + std_dev*randn( size(x) );